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yu  utilidor  height  (ft) 

6  effective  thickness  of  air  layer  (ft) 

tj0  overall  combined  plant  efficiency  (%) 

Axu  utilidor  wall  thickness  (ft) 
ua  kinematic  viscosity  of  air  (ftJ/s) 


tv 


LOSSES  FROM  THE  FORT  WAINWRIGHT 
HEAT  DISTRIBUTION  SYSTEM 


G.L.  Phetteplace,  W.  Willey  and  M.A.  Novick 


INTRODUCTION 

Fort  Wainwright,  located  near  Fairbanks,  Alaska,  is 
the  largest  U.S.  Army  base  located  in  an  extremely 
cold  climate.  Winter  temperatures  often  drop  to  -40° F, 
and  occasionally  temperatures  of  -50° F  or  lower  are 
experienced.  In  such  a  climate,  space  heating  can  no 
longer  be  considered  a  comfort;  it  becomes  a  true  ne¬ 
cessity.  The  supply  of  heat  must  be  continuous  and 
reliable.  Even  a  minor  outage  of  several  hours  can 
drop  indoor  temperatures  to  below  freezing. 

Fort  Wainwright,  like  many  large  military  bases,  has 
a  central  heat  and  power  plant.  This  plant  produces 
heat  and  electricity  to  meet  the  entire  needs  of  the 
base,  with  the  exception  of  some  outlying  buildings 
which  are  heated  by  other  means.  The  purpose  of 
this  report  is  to  examine  Fort  Wainwright’s  heat  dis¬ 
tribution  system  in  some  detail  in  order  to  estimate  its 
efficiency. 

At  Fort  Wainwright  during  the  winter,  heat  losses 
from  the  buried  heat  distribution  system  prevent  snow 
from  accumulating  over  it  in  all  but  the  coldest  periods. 
Figure  l  shows  the  bare,  dry,  ground  above  a  utilidor 
at  Fort  Wainwright  during  March  of  1979.  It's  not 
unusual  at  this  time  of  year  to  see  children  practicing 
baseball  for  the  coming  summer  months  on  bare  ground 
above  the  buried  utilidors.  After  witnessing  bare 
ground  over  nearly  all  the  sections  of  the  buried  heat 
distribution  system,  it's  hard  to  believe  that  the  heat 
losses  are  anything  but  staggering. 


This  apparent  waste  of  thermal  energy  prompted  a 
study  to  determine  its  general  magnitude.  By  using 
data  from  the  Fort  Wainwright  heat  and  power  plant, 
we  have  determined  the  magnitude  of  the  heat  loss  and 
have  developed  a  method  which  allows  us  to  rapidly 
evaluate  any  proposed  improvements  to  the  system. 

THE  FORT  WAINWRIGHT  HEAT 
AND  POWER  PLANT 

Before  going  into  detail  on  the  Fort  Wainwright 
heat  and  power  plant,  let’s  consider  some  of  the  gener¬ 
al  characteristics  of  combined  heat  and  power  plants. 

Conventional  production  of  electricity  requires  that 
fuel  energy  first  be  converted  into  heat  before  it  can 
be  converted  into  mechanical  energy  and  eventually 
into  electricity.  The  overall  efficiency  of  this  process 
of  changing  fuel  to  electricity  is  seldom  over  35%  in 
the  most  modern  plants,  and  frequently  less  than  30% 
in  the  older  facilities.  The  limiting  factor  is  the  con¬ 
version  of  heat  to  mechanical  energy,  in  which  only  a 
fraction  of  the  heat  may  be  used.  The  availability  of 
the  remaining  heat  is  too  low  and  it  must  be  rejected. 
In  a  conventional  power  plant,  this  waste  heat  is  re¬ 
jected  to  the  environment  by  some  means.  In  a  com¬ 
bined  heat  and  power  plant  the  waste  heat  is  used  to 
provide  heat  for  distribution  to  buildings.  The  waste 
heat,  which  is  a  liability  in  a  conventional  plant,  now 
becomes  an  asset  in  a  combined  heat  and  power  plant. 


Figure  7.  The  ground  above  a  buried  utllidor  in  March  at  Fort  Wainwright. 


Table  1.  Heat  balances  for  typical  electrical  and  com¬ 
bined  heat  and  power  plants. 


Output  us  percentage  ot 
fuel  energy  input 

t  lee  tricity- 
only  plant 

Combined  heat  and 
power  plant 

Electricity 

35 

25 

Waste  heat 

50 

0 

Useful  heat 

0 

60 

Stack  losses 

10 

10 

Parasitic  and  plant  heat  loss 

5 

5 

Total 

100 

100 

The  recovery  and  use  of  this  waste  heat  is,  however, 
not  without  penalty.  Normally,  the  heat  is  rejected 
at  a  temperature  too  low  for  most  space-heating  appli¬ 
cations.  When  the  heat  is  rejected  at  higher  tempera¬ 
tures,  where  it  can  be  used,  a  portion  of  electrical 
production  is  lost.  Overall,  however,  the  net  efficiency 
of  the  plant  is  greatly  increased.  Consider  the  heat  bal¬ 
ances  in  Table  1 ,  which  might  be  typical  of  each  type 
of  plant  during  the  heating  season. 

In  the  case  presented  in  Table  1  the  useful  output 
of  the  plant  increased  from  35  to  85%.  The  thermal 


energy  is  always  of  lesser  value,  however,  and  can  usu¬ 
ally  be  produced  for  about  25%  of  the  cost  of  an  equiv¬ 
alent  amount  of  electrical  energy.  Even  so,  the  total 
value  of  the  plant  output  relative  to  electrical  energy 
has  increased  from  35  to  40%  1 25  +  0.25(60) )  which 
represents  at  14%  increase.  During  non-heating  peri¬ 
ods  a  properly  designed  heat  and  power  plant  would 
be  able  to  deliver  the  same  efficiency  as  a  plant  that 
produces  electricity  only. 

As  illustrated  by  this  example,  the  economic  bene¬ 
fits  of  combined  heat  and  power  production  arc  modest. 


The  relative  energy  savings,  however,  are  of  much 
greater  magnitude.  If  we  assume  that  the  combined 
heat  and  power  plant  supplies  all  the  heat  required, 
we  can  determine  the  total  fuel  requirements  for  non 
central  heat  production  as  follows: 

%  heat  required  =  %  of  fuel  input  for  heating 
from  a  combined  plant 

-  %  of  energy  lost  in  distrubution. 

Assuming  that  10%  of  the  heat  is  lost  in  the  heat  dis¬ 
tribution  system  and  using  the  60%  heat  output  value 
for  a  combined  plant  from  Table  1 ,  we  have 

Heat  required  =  60-0.10(60)  =  60-6  =  54%. 

In  a  non-central  heating  scheme  this  would  be  supplied 


the  turbines.  The  plant  has  five  turbine-generator  sets 
with  a  total  nameplate  capacity  of  22  MW  electricity. 
Three  of  the  turbines  are  5-MW  extraction-condensing, 
one  is  a  5-MW  extraction-backpressure,  and  the  re¬ 
maining  one  is  a  2-MW  extraction-condensing  type. 
Basically,  in  an  extraction  turbine  the  steam  is  expand¬ 
ed  partially  to  a  lower  pressure  where  a  fraction  is  ex¬ 
tracted  from  the  turbine  for  space  heating  or  industrial 
processes.  The  remaining  steam  continues  through  the 
low-pressure  stages  of  the  turbine.  If  the  steam  leaves 
the  turbine  below  atmospheric  pressure  and  is  subse¬ 
quently  condensed,  the  turbine  is  called  an  extraction- 
condensing  turbine.  If  the  steam  leaves  the  turbine  at 
higher  pressures  and  is  used  for  other  purposes  the  tur¬ 
bine  is  called  an  extraction-backpressure  turbine.  Con¬ 
ventional  (non-extraction)  turbines  may  be  of  either 
the  backpressure  or  condensing  type.  At  Fort  Wain- 
wright  the  extraction  occurs  at  approximately  100 


by  individual  heating  plants.  Their  efficiency  is  some¬ 
what  dependent  on  fuel.  Assuming  oil  is  the  fuel,  70% 
is  a  typical  efficiency.  Thus,  the  fuel  energy  required 
for  the  individual  heating  plant  would  be 

Fuel  required  for  individual  heating  plants 

=  54/0.70  =  77%  of  central  plant. 

Electric  generation  is  somewhat  more  efficient  for 
the  electric-only  plants.  It  would  require  only  71% 

1 1 00%(25%/35%)  1  of  the  fuel  required  by  the  com¬ 
bined  heat  and  power  plant  to  generate  the  same 
amount  of  electricity.  Thus  the  overall  fuel  use  for 
non-central  heating  (77%)  and  central  electricity  pro¬ 
duction  is  1.48  times  that  of  combined  heat  and  power 
production.  Thus  the  energy  saving  is  nearly  50%. 
Central  heating  also  has  the  added  advantage  of  being 
able  to  use  less  expensive  or  more  plentiful  fuels  (i.e. 
coal,  heavy  oil,  nuclear,  as  opposed  to  fuel  oil  and 
natural  gas). 

Now  that  we  have  discussed  a  few  of  the  advantages 
of  combined  heat  and  power  plants,  let's  discuss  the 
Fort  Wainwright  plant  in  some  detail.  Figure  2  is  a 
simplified  flow  diagram  for  the  plant  (from  Rossie  et 
al.  1975). 

Steam  is  generated  in  six  of  the  eight  boilers,  as  the 
two  original  boilers  of  the  plant  are  not  operational. 
Each  of  the  operating  boilers  is  rated  at  1 50,000  Ibm 
of  s team  per  hour.  The  steam  leaves  the  boilers  at  a 
pressure  of  400  to  420  psig,  is  superheated  to  a  tem¬ 
perature  between  650  and  7S0°F,  and  is  then  fed  to 


psig  in  all  the  turbines.  The  backpressure  turbine  ex¬ 
hausts  at  10  psig  and  this  steam  is  used  for  power  plant 
heating. 

The  overall  combined  efficiency  of  both  heat  and 
electrical  generation  at  Fort  Wainwright  can  be  calcu¬ 
lated  based  on  records  maintained  at  the  plant.  The 
overall  combined  efficiency  is  simply 

r?0  =  Overall  combined  plant  efficiency 

_  Steam  heat  out  and  electricity  out 
fuel  energy  input 

Consider  the  following  sample  calculation  for  the 
month  of  January  1975.  From  steam  tables  we  can 
findthe  enthalpy  of  350°F  steam  at  105  psia  ('  90 
psig)  to  be  1204.5  Btu/lbm.  Similarly,  we  find  the 
enthalpy  of  the  condensate  to  be  118  Btu/lbm.  The 
difference  between  these  two  values  represents  the 
heating  energy  sent  out  per  pound-mass  of  steam, 

1086.5  Btu/lbm.  For  January  1975,  1.33x10®  Ibm 
of  steam  was  sent  out,  along  with  9.05x1 06  kWh  of 
electricity.  In  the  same  period  27,604  tons  of  coal 
were  used.  If  the  average  heat  content  of  the  coal  is 
assumed  to  be  8600  Btu/lbm,  then  the  efficiency  r?0 
is  calculated  as  37.0%.* 

Notice  that  this  is  less  than  half  the  efficiency  of 
the  hypothetical  plant  discussed  earlier.  Normally  the 
Fort  Wainwright  heat  and  power  plant  runs  at  a  com¬ 
bined  efficiency  of  about  60%.  The  values  are  plotted 
for  the  combined  efficiency  over  a  four-year  period 
in  Figure  3.  Even  a  60%  efficiency  is  significantly  lower 
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Figure  3.  Overall  combined  plant  efficiency  and  normalized  electric  output 
for  the  Fort  Wainwright  power  plant. 


Figure  4.  Steam  flow  and  degree-day  data  for  the  Fort  Wainwright  power 
plant. 


than  possible,  primarily  due  to  the  design  of  the  plant 
itself.  Steam  pressure  is  too  low  for  efficient  electric 
generation,  and  the  plant  was  also  designed  for  high 
reliability  which  causes  some  compromise  of  efficiency. 
It  still  supplies  heat  and  electricity  more  economically 
and  reliably  than  any  other  method. 

As  can  be  seen  from  Figure  3,  the  electrical  load  of 
the  power  plant  peaks  in  the  winter  months,  as  would 
certainly  be  expected  for  any  electrical  utility  in  such 
a  climate.  Figure  4  illustrates  the  steam  flow  from  the 
plant  over  the  same  four-year  period.  Also  shown  in 
Figure  4  arc  heating  degree-day  data  for  this  period. 

The  degree-day  is  a  measure  of  heating  requirements, 
and  as  suspected,  the  steam  flow  from  the  plant  fol¬ 
lows  the  degree-day  data  closely.  The  one  discrepancy 
occurs  in  the  summer  months,  when  heat  demand  is 


low  or  nonexistent.  Part  of  this  discrepancy  is  due  to 
the  fact  that  domestic  hot  water  is  heated  with  steam 
from  the  distribution  system  and  thus  the  load  from  it 
is  relatively  constant.  Based  on  the  number  of  residents 
and  employees  on  the  base,  domestic  hot  water  heating 
could  certainly  not  account  for  more  than  5%  of  the 
maximum  system  load.  The  remaining  gap  of  about 
20%  is  mostly  due  to  heat  losses  from  the  piping  sys¬ 
tem.  These  will  be  examined  in  detail  below. 


HEAT  DISTRIBUTION  AT 
FORT  WAINWRIGHT 

Steam  extracted  from  the  turbines  at  100  psig 
(actually,  the  system  has  been  using  90-psig  steam 


5 


recently)  is  distributed  to  nearly  all  the  buildings  on 
base  to  meet  both  space  heating  and  domestic  hot 
water  heating  needs.  In  85  to  90%  of  the  buildings, 
the  steam  is  used  directly  for  space  heating  after  pass¬ 
ing  through  a  pressure  reducing  valve.  The  remaining 
buildings  have  hot  water  heating  systems  and  a  water  - 
to-steam  heat  exchanger. 

The  distribution  system  itself  consists  primarily  of 
insulated  steel  pipes  in  utilidors.  Several  exceptions 
are  pipes  which  are  directly  buried  in  the  ground.  The 
utilidors  range  in  size  from  1  x  1  ft  (inside  dimensions) 
up  to  7x  9  ft.  The  larger  sizes  can  be  easily  walked  in 
for  pipe  servicing.  The  total  length  of  the  distribution 
system  is  about  28  miles  with  a  burial  depth  of  be¬ 
tween  2  and  6  ft. 

Nearly  all  of  the  piping  is.insulated,  but  much  of 
the  insulation  was  damaged  during  the  flood  of  Aug¬ 
ust  1967.  Part  of  the  damaged  insulation  was  re¬ 
placed,  and  in  other  instances,  insulation  was  added 
to  the  original  insulation  to  restore  its  thermal  resis¬ 
tance  to  the  original  value. 

Complete  maps  of  the  Fort  Wainwright  heat  dis¬ 
tribution  system  were  updated  in  1977  and  are  there¬ 
fore  quite  accurate.  Using  these  maps,  we  determined 
the  length  and  size  of  each  pipe  segment  and  tabulated 
this  information  under  the  appropriate  category  of 
utilidor  and  pipe  sizes.  From  this  table  the  total 
lengths  of  each  configuration  of  utilidor,  supply  pipe 
and  return  pipe  size  were  determined.  Appendix  A 
contains  a  table  of  some  200  different  combinations. 


HEAT  LOSSES  FROM 
UTILIDOR  SYSTEMS 

The  analysis  of  heat  losses  from  utilidors  is  not  a 
simple  problem.  Heat  flows  from  the  warm  pipes 
through  a  series  of  thermal  resistances  to  the  air  at 
the  ground  surface.  The  resistances  involved  are: 

1 .  Convective  heat  transfer  from  the  flow¬ 
ing  fluid  to  the  pipe’s  inner  surface. 

2.  Conductive  heat  transfer  through  the 
pipe  wall  to  the  outer  pipe  surface. 

3.  Conductive  heat  transfer  through  the 
insulation  to  its  outer  surface. 

4.  Combined  convective,  conductive  and 
radiative  heat  transfer  from  the  insulation  sur¬ 
face  to  the  inner  wall  of  the  utilidor. 

5.  Conductive  heat  transfer  through  the 
utilidor  wall. 

6.  Conductive  heat  transfer  through  the 
soil  to  the  ground  surface. 

7.  Finally,  convective  heat  transfer  to  the 
air  at  the  ground  surface. 


Fortunately,  resistances  1,  2  and  7  are  very  small 
in  comparison  to  the  others  and  can  be  neglected  with 
little  error.  Unfortunately,  resistance  4  is  very  compli¬ 
cated.  It  will  be  treated  in  a  simplified  manner. 

Consider  the  most  complicated  case,  both  a  steam 
and  condensate  (return)  pipe  inside  the  same  utilidor. 
(Actually,  instances  do  occur  where  several  steam 
and/or  condensate  pipes  may  share  a  common  utilidor. 
In  these  instances,  the  supply  and  return  pipes  are 
treated  as  pairs  and  additional  pipes  are  treated  as  single 
pipes.  The  total  heat  loss  is  then  assumed  to  be  the 
sum  of  all  the  pipes  included.)  The  overall  heat  flow 
from  the  warm  pipes  is  given  as 


where  q  =  overall  heat  loss  per  unit  length  (Btu/hr  ft) 

Q  =  overall  subsystem  heat  loss  (Btu/hr) 

L  =  subsvstem  length  (ft) 

U0  =  overall  t/-value  of  subsystem  (Btu/hr  ft  °) 

Ts  =  steam  temperature  (°F) 

7g  =  ground  surface  temperature  (°F). 

The  overall  conductance  value,  U0  is  simply  : 

^o  =  1/«o  (2) 

where  R0  is  the  overall  thermal  resistance  (hr  ft  °F/ 
Btu). 

Since  the  heat  flow  is  basically  a  series  phenomenon, 
we  can  simply  add  the  individual  resistances  to  heat 
flow  to  obtain  the  overall  resistance  to  heat  flow  R0 : 

*o  =  *ei  +  Kea  +  *u  +  Ks  (3) 

where  /?ei  =  effective  thermal  resistance  of  the  in¬ 
sulation  on  the  pipes  (hr  ft  °F/Btu) 

Rea  =  effective  thermal  resistance  of  the  air 
within  the  utilidor  (hr  ft  °F/Btu) 

Ru  =  thermal  resistance  of  the  utilidor  (hr  ft 
°F/Btu) 

R s  =  thermal  resistance  of  the  surrounding  soil 
(hr  ft  °F/Btu) 

To  determine  the  effective  thermal  resistance  of  the 
pipe  insulation  Rej  we  must  consider  the  parallel  heat 
flow  from  each  of  the  pipes  through  its  respective  in¬ 
sulation.  The  sum  of  the  two  heat  flows  will  be  the 
total  heat  flow  q: 


yr, i  t  rc-7-ci 


*si 


(4) 
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where  7V  =  outer  surface  temperature  of  the  steam 
pipe  insulation  (°F) 

7cj  =  outer  surface  temperature  of  the 
condensate  pipe  insulation  (°F) 

R%i  =  thermal  resistance  of  the  steam  pipe  in¬ 
sulation  (hr  ft  °F/Btu) 

/?ci  =  thermal  resistance  of  the  condensate  pipe 
insulation  (hr  ft  °F/Btu) 
rc  =  condensate  return  temperature  (°F). 


and  similarly  for  the  condensate  pipe  insulation 

rmi  =  (7c+rci)/2-  01) 

Next  we  need  to  find  the  effective  resistance  of  the 
air  Rel.  First,  consider  the  heat  flow  equation  for  the 
air  layer: 

Q=l'lRjTirTu.)HMRjTcrTui)  (12) 


The  total  heat  flow  can  also  be  written  as 


(5) 


By  combining  eq  4  and  5  we  can  find  the  effective 
thermal  resistance  of  the  pipe  insulation  Rtl  as 


*ei  = 


Tci 


*si 


V  ^  *ci 


(6) 


where  RSi  =  thermal  resistance  between  the  steam 

pipe  insulation  surface  and  the  inner  wall 
of  the  utilidor  (hr  ft  °F/Btu) 

Rci  =  thermal  resistance  between  the  condensate 
pipe  insulation  surface  and  the  inner  wall 
of  the  utilidor  (hr  ft  °F/Btu) 
ruj  =  temperature  of  the  inner  wall  of  the  utilidor 

(°F). 

The  heat  flow  can  also  be  expressed  in  terms  of  the 
effective  thermal  resistance  of  the  air  /?ea  as 


and  the  individual  pipe  insulation  resistance  can  be  q-  1//?ea  (7V  -  ruj).  (1 3) 

calculated  from  the  following  expression  (Holman 

1972)  By  combining  eq  12  and  13  we  find 


ln[(Ds+2Xsi)/Ps) 

2nkx 


(7) 


and 


*ci 


ln[(Oc+2Xci)/Pc] 

2nk- 


(8) 


where  Ds  =  outside  diameter  of  the  steam  pipe  (ft) 
Xsj  =  insulation  thickness  on  steam  pipe  (ft) 

6 1  =  insulation  thermal  conductivity  (Btu/hr 
ft  °F) 

Oc  =  outside  diameter  of  the  condensate  pipe 
(ft) 

2fcj  =  insulation  thickness  on  the  condensate 
Pipe  (ft). 


The  thermal  conductivity  of  insulation  is  normally 
a  function  of  its  mean  temperature.  For  calcium  sili¬ 
cate  insulation,  using  data  from  Crocker  and  King  (1967), 
this  function  is  closely  approximated  by 


ki  =0.0221  +4.13x1(T5  Tmi  (9) 


where  7"mi  is  the  average  temperature  of  the  pipe  insula¬ 
tion  (°F).  The  temperature  of  the  steam  pipe  insulation  is 

^i  =  (Ts-rsi)/2  (10) 


1 _ _ 

i  ,  ?  lTcrrj 

^sa  ^ca  (f"si  “  ^ui) 


(14) 


Now  we  need  to  define  the  individual  thermal  resis¬ 
tances  caused  by  the  air  layer  £sa  and  RCi.  As  men¬ 
tioned  earlier,  the  heat  transfer  process  within  the  air 
space  is  coupled  convective,  conductive  and  radiative. 

A  full  treatment  is  beyond  the  scope  and  requirements 
of  this  study.  Instead,  we  will  treat  the  air  layer  as  a 
conductive  medium  and  determine  the  effective  ther¬ 
mal  conductivity  value  6ea  for  the  medium  based  on 
published  correlations  for  heat  transfer  in  concentric 
annuli. 

In  order  to  consider  this  region  as  an  annulus  we 
must  first  define  an  effective  diameter  for  the  rectangu¬ 
lar  utilidor.  The  effective  diameter  0eu  is  the  diameter 
of  a  circular  utilidor  which  would  have  the  same  inside 
surface  area  as  the  actual  rectangular  utilidor: 


^eu  =  ^  (*u  + 


(15) 


where  Xu  is  the  utilidor  width  (ft)  and  /u  the  utilidor 
height  (ft). 

The  effective  thermal  conductivity  for  the  air  within 
such  a  region  is  approximated  by  Gruber  ct  al.  (1961 ) : 

^a  =  0.40(/VG/Vp)°-20-^ 
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(16) 


r 


within  the  range  106  <  /VG/Vp  <  109 

where  kei  =  effective  thermal  conductivity  of  the  air 
(Btu/hr  ft  °F) 

Ng  =  Grashof  number  (dimensionless) 

Nf  =  Prandtl  number  (dimensionless) 

ki  =  thermal  conductivity  of  air  (Btu/hr  ft  °F). 

The  thermal  conductivity  of  air  k3  is  a  function  of 
its  temperature  and  can  be  approximated  as 

*a  =  0.01319  +  2.5x10‘s  fa  (17) 

where  Ta  is  the  bulk  air  temperature  within  the  an¬ 
nulus  (°F). 

The  Prandtl  number  of  air  Nf  is  also  a  function  of 
the  air  temperature  and  can  be  approximated  as 


to  evaluate  the  effective  thermal  conductivity  of  the 
air  kei,  we  can  find  the  thermal  resistance  due  to  the 
airspace,  it  is  simply 


■nlPeu/(P«  +  2*t|,M 

2  IT*., 


(24) 


for  the  steam  pipe.  For  the  condensate  pipe 


ln|Peu/(Pc  +  2Xcl)| 
2nkei 


(25) 


The  next  thermal  resistance  which  we  need  to  define 
is  that  of  the  utilidor  itself.  It  can  be  written  as 


/?„  = 


A*,, 


2  (*U  +  Yu)ku 


(26) 


Np  =0.71 85  -1.275x1  (T4  Tr 


(18) 


The  bulk  air  temperature  T3  can  be  approximated 
by  the  average  of  the  effective  insulation  surface  tem¬ 
perature  Tt  j  and  the  utilidor's  inner  surface  tempera- 


ture7ui: 


ra  =  (rei  +  7iii)/2- 


(19) 


The  effective  insulation  surface  temperature  is 
given  as  the  weighted  average  of  the  steam  and  con¬ 
densate  insulation  surface  temperatures: 


^i  = 


W+2*,i)+7ci(Pc+2Xc 

Ds+2Xsi+Dc  +  2X£j 


(20) 


The  Grashof  number  NQ  given  in  eq  16  is  repre¬ 
sented  for  this  case  as 


"g  = 


g^ei-^K*3) 

(ra+459.7)(ua)2 


(21) 


where g  =  gravitational  constant  (ft/s2) 

6  =  effective  thickness  of  air  layer  (ft) 
ua  =  kinematic  viscosity  of  air  (ft2/s). 

The  effective  thickness  of  the  air  layer  6  is  given  by 

0eu-(0s+2*si  +  Dc  +  2Xci) 


6  = 


(22) 


The  kinematic  viscosity  of  air  ua  is  also  a  function 
of  the  air  temperature  Ta.  It  can  be  approximated  as 

v3  =  1 .26x10"*  +  5.4x10" 7  ra.  (23) 


where  AXU  is  the  utilidor  wall  thickness  (ft)  and  ku  is 
the  thermal  conductivity  of  the  utilidor  material  (Btu/ 
hr  ft  °F). 

And,  finally,  the  last  thermal  resistance  we  need  to 
define  is  that  of  the  soil  system: 


1/*sS 


(27) 


where  kf  is  the  thermal  conductivity  of  the  soil  (Btu / 
hr  ft  °F)  and  S  is  the  shape  factor  of  the  utilidor  (di¬ 
mensionless). 

For  a  rectangular  utilidor  the  shape  factor  S  is  given 
(from  Holman  1972)  as 


5=HM'4:r4:r' 


v -0.078, 


(28) 


where  is  the  utilidor  burial  depth  (ft). 

Now  we  have  expressions  for  all  the  thermal  resis¬ 
tances  given  in  eq  3.  With  these  we  can  find  the  heat 
loss  from  the  utilidor.  The  calculation  procedure  is 
not  completely  straightforward,  however. 

Notice  that  both  the  effective  thermal  resistance  of 
the  pipe  insulation  Rci  and  the  effective  thermal  resis¬ 
tance  of  the  air  space  Rt3  arc  functions  of  tempera¬ 
tures  which  arc  not  initially  known.  It  is  necessary, 
therefore,  for  us  to  solve  this  problem  iteratively.  To 
do  so  we  first  guess  at  the  unknown  required  tempera¬ 
tures.  We  then  use  the  procedure  outlined  to  calculate 
the  heat  flow  q.  Given  this  result  we  can  calculate  the 
outer  surface  temperature  of  the  utilidor  7"uo  by 

7uo=  V9*s  <29> 


Now  that  wc  have  all  of  the  necessary  expressions 


and  similarly  the  inner  surface  temperature  of  the 
utilidor  is  given  as 


START 


HEAD  IH 

NUMBER  0* 
SECTIONS  <-Nl 


READ  IN  SET  Of 
UTILUJOH  ANO 
NH  DIMENSIONS 


UNO  MEAT  VOSS  AMD 
RETURN  INSL  SURFACE 
TIME  BASED  ON 
CONDENSATE  T£MP 


TlNO  HEAT  LOSS  AND 

Surely  inSv  sure  ac  i 

I  TEMP  tASCO  ON 
STEAM  TEMP 


1 1*0  h«  at  vosseAtto 
ON  ST l  AM  T*MP 
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RESISTANCE  VALUES 


'  RETURN  ^ 
.PIPE  D'A-  »  0 


r  SUPPLY  PIPE  ONLY 
I  CASE  •  I 


f  <  MO  SU*i  r  A  AO  *  t  TURN 
INSt  SU ft  I  ACE  TEAMS 
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AT  THE  IN  Si  SURE  ACC 


/  SUPPLY  ^ 
v/ME  OfA  «  Q  . 


SUPPLY  PIPE  0|A  »  0 


/  ABSOLUTE  N, 

'error  in  neat  loss 

\  .  o  tois  / 


[ElNO  HEAT  LOSS 

Over  entire 
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Figure  5.  Computer  program  flow  chart. 


INSL  SURFACE  TfMB 
EQUALS  THAT  OF 
SURAL Y  RIPE  INSV. 


insl.  surface  temp, 
equals  that  of 

RETURN  PIPE  INSL. 


EFFECTIVE  INSI  SURFACE  I 
TEMP  EQUALS  WEKJMTfO  j 
AVC  Of  SUPPLY  ANO  j 
RETURN  PIPE  INSL  TEMPS.  1 
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(30) 


7"ui  ^uo+<?^U' 

Now  we  need  to  recalculate  the  insulation  surface 
temperatures  7"SI  and  Tci.  We  can  do  this  by  consider¬ 
ing  a  heat  balance  at  the  insulation  surface  and  using 
the  new  value  for  the  utilidor  inner  surface  temperature 
Tu |  just  calculated.  The  heat  balance  gives  the  follow¬ 
ing  equations: 

j  - - — - + - — -  (31) 

Sl  1  +  («si/«sa)  1  +  (KSa/*si)  ' 

and 

t  =  _ — h - + - -  (32) 

Ci  l+(*ci/«ca)  1  +  (*ca/*ci)‘ 

Now  we  can  use  these  new  temperatures  to  recal¬ 
culate  the  heat  loss  q.  If  the  heat  loss  is  close  to  the 
value  previously  calculated,  we  have  found  the  answer; 
if  it  is  not  we  must  recalculate  the  temperatures  and 
heat  loss  again.  This  process  continues  until  the  heat 
loss  value  has  stabilized  and  the  answer  is  found.  A 
similar,  but  simpler,  procedure  can  be  used  for  single 
pipes  in  utilidors. 

A  computer  program  written  to  accomplish  this 
calculation  procedure  for  an  entire  heat  distribution 
system  is  listed  in  Appendix  B.  A  flow  chart  of  the' 
program  is  given  in  Figure  5.  In  Appendix  C,  sample 


input  and  output  data  for  the  program  are  given. 

The  advantages  of  computer  implementation  of  this 
calculation  scheme  are  obvious.  As  well  as  allowing  us 
to  determine  the  heat  loss  of  any  heat  distribution  sys¬ 
tem  using  utilidors,  it  allows  us  to  evaluate  modifica¬ 
tion  to  the  system  or  operating  parameters  to  deter¬ 
mine  their  effect.  In  the  next  section  we  will  examine 
the  results  for  the  Fort  Wainwright  heat  distribution 
system.  Modifications  and  their  effect  will  also  be 
studied. 


SYSTEM  EFFICIENCY  AND 
POSSIBLE  IMPROVEMENTS 

Using  a  computer  program  following  the  method 
outlined  above,  we  have  found  the  total  heat  loss  from 
the  distribution  system  to  be  2.045x10s  MBtu/yr. 

This  value  assumes  that  the  average  annual  air  tempera¬ 
ture  of  25. 7°  F  can  be  used  as  the  ground  surface  tem¬ 
perature.  It  also  assumes  that  all  pipes  are  insulated 
with  1  in.  of  calcium  silicate  insulation. 

The  heat  loss  from  the  utilidor  system  actually  var¬ 
ies  over  the  course  of  the  year  due  to  the  fluctuations 
in  the  ground  surface  temperature.  Figure  6  shows 
how  the  air  temperature  varies  over  an  average  year 
as  well  as  its  effect  on  the  heat  loss.  The  temperatures 
used  in  this  case  represent  the  average  monthly  air  tem¬ 
peratures  for  the  Fairbanks  area. 
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Figure  6.  Ground  surface  temperature  and  heat  losses  for  an  average 
year. 


For  a  more  specific  case,  calendar  year  1978,  let’s 
compare  the  heat  loss  over  the  course  of  the  year  with 
the  heat  supplied  to  the  distribution  system.  We  can 
express  the  heat  supplied  to  the  system  over  the  course 
of  the  year  in  terms  of  British  Thermal  Units  by  recall¬ 
ing  that  each  pound-mass  of  steam  represents  1086.5 
Btu  of  heat  at  the  average  steam  and  condensate  condi¬ 
tions.  Table  2  summarizes  the  plant  input  to  the  sys¬ 
tem  as  well  as  the  calculated  losses. 

The  estimated  percentage  heat  loss  during  the 


month  of  August  is  more  than  double  that  of  the 
month  of  January,  based  in  each  case  on  the  heat  ex¬ 
ported  from  the  power  plant.  This  is  due  to  the  small 
heat  loss  fluctuations  over  the  yearly  cycle  compared  to 
the  steam  flow  from  the  plant  that  varies  to  a  greater 
extent.  If  we  normalize  each  of  these,  that  is  divide 
them  by  the  largest  value  over  the  year,  this  effect  is 
clearly  illustrated.  Figure  7  gives  the  normalized  heat 
export  by  the  plant  as  well  as  the  normalized  heat  loss 
for  the  1978  calendar  year. 


Figure  7.  Normalized  heat  export  and  heat  loss  for  1 9  78. 


Table  2.  Exported  heat  and  calculated  heat  losses  for  calendar  year  1978. 


Rate  of  heat  loss 

Normalized 

Normalized 

Steam 

A  verage  rate 

A  verage  air 

Calculated  rate 

as  percentage  of 

rate  of  heat 

rate  of 

exported 

of  heat  export 

temperature 

of  heat  loss 

rate  of  heat  exported 

export 

heat  loss 

Month 

Days 

(hlbrn) 

(MBtu/hr) 

(°F) 

(MBtu/hr) 

(%) 

(%) 

(%) 

|an 

31 

151,633 

221.4 

0.1 

25.4 

11.5 

100.0 

100.0 

Feb 

28 

125,540 

203.0 

3.9 

25.1 

12.4 

91.7 

98.8 

Mar 

3t 

127,576 

186.3 

14.0 

24.3 

13.0 

84.1 

95.7 

Apr 

30 

100,559 

151.7 

34.8 

22.6 

14.9 

68.5 

89.0 

May 

31 

73,413 

107.2 

50.2 

21.3 

19.9 

48.4 

83.9 

)une 

30 

63,516 

95.8 

S4.6 

21.0 

21.9 

43.3 

82.7 

)uly 

31 

59,134 

86.4 

63.5 

20.2 

23.4 

39.0 

79.5 

Aug 

31 

57,309 

83.7 

59.5 

20.6 

24.6 

37.8 

81.1 

Sept 

30 

66,324 

100.1 

46.8 

21.6 

21.6 

45.2 

85.0 

Oct 

31 

107,812 

157.4 

23.3 

23.5 

14.9 

71.1 

92.5 

Nov 

30 

125,039 

188.7 

8.6 

24.7 

13.1 

85.2 

97.2 

Oec 

31 

137,662 

201.0 

3.3 

25.2 

12.5 

90.8 

99.2 

Total 

365 

1,195,517 

- 

- 

~ 

- 

- 

Averages 

30.4 

99,626 

148.6 

30.2 

23.0 

15.5 

67.1 

90.4 
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Table  3.  The  effect  of  increasing  insulation  thickness  on  all  pipes. 


Steam  pipe 

Condensate  pipe 

Total 

Change  in 

Total  ins. 

Change  in 

Heat  savings  per  unit 

ins.  thickness 

ins.  thickness 

heat  loss 

heat  loss 

volume 

ins.  volume 

of  insulation  added 

Cose 

(in.) 

(in.) 

(MBtu/yr) 

(MBtu/yr) 

(ft*) 

(ft*) 

(MBtu/yr  ft') 

Reference 

1 

1 

2.045  X 1 0! 

1.202x10s 

1 

1 

2 

1.920X105 

-  1.25x10* 

1. 788x10s 

5.860X10* 

0.2133 

2 

2 

1 

1.551x10s 

-4.940x10* 

1.817X10* 

6.150X10* 

0.8033 

3 

2 

2 

1.430X10’ 

-6.150X10* 

2.403X10* 

1.201x10s 

0.5121 

4 

3 

1 

1.317x10s 

-7.280X10* 

2.432X10* 

1.230x10s 

0.5919 

5 

3 

2 

1.198x10s 

-8.470X10* 

3.019X10* 

1.817x10s 

0.4662 

Table  4.  The  effect  of  increasing  insulation  thickness  of  pipes  in  utilidors  3Vi  ft  by  3’/2  ft  or  larger. 

Steam  pipe 

Condensate  pipe 

Total 

Change  in 

Total  ins. 

Change  in 

Heat  savings  per  unit 

ins.  thickness 

ins.  thickness 

heat  loss 

heat  loss 

volume 

ins.  volume 

of  insulation  added 

fuse  (in.) 

(in.) 

(MBtu/yr) 

(MBtu/yr) 

<ft') 

<ft'> 

( MBtu/yr  ft') 

Reference 

i 

1 

2.045X10* 

- 

1.202X10* 

- 

- 

l 

i 

2 

1.940X1 0s 

-1.050X10* 

1.649X1 0s 

4.470X10* 

0.2349 

2 

2 

1 

1.621X10* 

-4.240X10* 

1.677X10* 

4.750X10* 

0.8926 

3 

2 

2 

1.518x10s 

-5.270X10* 

2.125X10* 

9.230X10* 

0.5710 

4 

3 

1 

1.419X10* 

-6.260X10* 

2.152X10* 

9.500X10* 

0.6589 

5 

3 

2 

1.317X10* 

-7.280X10* 

2.600X10* 

1.398X10* 

0.5207 

As  stated  earlier,  the  heat  loss  for  an  average  year 
is  approximately  2.045x  10s  MBtu.  During  1979  the 
cost  of  heat  at  Fort  Wainwrjght  was  about  $4.60/ 
MBtu.  Thus,  the  yearly  cost  of  heat  loss  is  about 
$940,000.  With  energy  costs  rising  as  they  are,  this 
cost  will  soon  be  over  one  million  dollars  per  year. 

Let  us  now  consider  some  possible  improvements  to 
the  system  and  their  effect  on  heat  losses. 

The  most  obvious  method  of  reducing  heat  losses 
is  to  increase  insulation  thickness  on  the  pipes.  This 
is,  however,  a  very  expensive  proposition.  Although 
we  will  not  try  to  assess  the  cost  of  increasing  insula¬ 
tion  thicknesses,  we  will  provide  data  on  its  anticipated 
effect  on  heat  losses  from  the  system.  Table  3  gives 
the  results  of  several  sets  of  calculations  with  varying 
insulation  thicknesses.  Also  given  are  the  volumes 
and  the  increases  over  the  current  amount. 

In  the  last  column  of  Table  3  the  heat  savings  per 
unit  of  insulation  added  are  given.  These  quantities 
in  conjunction  with  the  cost  of  insulation  would  allow 
one  to  examine  the  relative  economics  of  each  case. 

As  can  be  seen  from  the  table,  case  2  offers  the  best 
opportunity  for  being  economically  viable.  The  heat 
savings  per  unit  of  insulation  added  arc  also  relatively 
high  for  case  4.  In  each  of  these  cases,  all  additional 
pipe  insulation  was  added  to  the  steam  pipes  alone. 

Another  possibility  investigated  was  insulating  only 
the  pipes  in  the  larger  size  util  idors.  Access  to  these 


pipes  would  be  easier  since  these  utilidors  can  be  en¬ 
tered,  as  opposed  to  smaller  utilidors  which  would  need 
to  be  excavated;  thus  the  economics  for  this  situation 
should  be  more  attractive.  Table  4  gives  results  similar 
to  Table  3  except  that  only  pipes  in  utilidors  3'/2X  3 Vi  ft 
or  larger  received  additional  pipe  insulation. 

As  before,  case  2  is  still  the  most  attractive  and  case 
4  is  the  next  best.  Also  notice  that  in  each  case  the 
heat  saving  per  unit  of  insulation  added  is  greater  than 
for  the  same  case  when  all  the  pipes  were  reinsulated. 
This,  coupled  with  the  fact  that  these  pipes  in  the  larg¬ 
er  utilidors  should  be  easier  to  insulate,  indicates  case  2 
as  the  best  approach.  In  addition,  it  seems  likely  that 
larger  pipes,  which  are  normally  in  the  larger  utilidors, 
would  be  easier  to  insulate  in  themselves.  This  would 
result  in  still  lower  insulating  costs. 

Another  method  which  could  be  used  to  reduce  heat 
losses  is  to  reduce  the  distribution  temperature  of 
steam.  Currently  steam  conditions  are  approximately 
90  psig  and  375°F,  although  the  saturation  tempera¬ 
ture  of  90  psig  steam  is  about  331°F.  Thus,  the 
steam  has  about  44°F  of  superheat.  If  this  amount  of 
superheat  and/or  the  distribution  pressure  were  re¬ 
duced  (i.e.  saturation  pressure  of  steam  reduced),  the 
heat  losses  could  also  be  reduced.  This  option  was 
studied  only  to  the  extent  of  determining  what  the 
resulting  heat  losses  would  be  for  any  given  steam  tem¬ 
perature.  Figure  8  shows  the  effect  of  steam  tempera¬ 
ture  on  the  total  system  heat  losses. 
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Steam  Temperature  (°F) 

Figure  8.  System  heat  losses  as  a  function  of  steam  temperature. 


We  feel  that  insulating  the  steam  pipes  in  the  larger 
utilidors  and  reducing  the  steam  temperature  hold  the 
most  promise  for  reducing  heat  losses  from  the  Fort 
Wainwright  system.  Although  other  possibilities  may 
exist,  they  will  not  be  studied  here.  These  two  methods 
would  produce  the  following  cost  savings. 

1.  Additional  pipe  insulation.  Based  on  the  limited 
input  data,  the  most  promising  alternative  would  be  to 
add  insulation,  probably  on  the  order  of  an  inch  or  so, 
to  the  steam  pipes  only.  It  is  most  likely  that  this  al¬ 
ternative  would  be  the  most  attractive  from  an  econom¬ 
ic  viewpoint.  Using  10%  as  the  time  value  of  money 
and  5%  as  the  escalation  rate  for  coal  (see  U.S.  Army 
Corps  of  Engineers  1978)  we  can  calculate  the  present 
worth  of  the  heat  savings  resulting  from  the  addition 
of  insulation.  For  case  2  from  Table  3  (addition  of 
one  inch  of  insulation  to  the  steam  line  only)  the  re¬ 
duction  in  heat  loss  is  4.94x10*  MBtu/yr.  For  1979 
the  heat  cost  at  Fort  Wainwright  was  $4.60/MBtu 
(from  Flanders  and  Coutts,  in  prep.).  Of  this  $2.10  per 
MBtu  is  attributable  to  fuel  costs  with  the  remaining 
$2.50  per  MBtu  attributable  to  ownership,  operations 
and  maintenance  costs.  If  we  assume  that  the  insula¬ 
tion  added  will  have  a  useful  lifetime  of  20  years,  the 
present  value  of  the  heat  savings  is  $2,307,770*  where 
the  multiplication  factors  (12.11  and  8.51)  are  the 


*  ”  -  4'94*’°4  MBtuJ(2.IO  ^  X  12.1 t)  ♦  (2.50  ^  x  8.51 


present  worth  factors  for  escalating  fuel  cost  savings 
and  fixed  plant  savings,  respectively. 

This  represents  the  maximum  amount  which  could 
be  justified  for  adding  one  inch  of  insulation  to  the 
steam  pipes  only. 

2.  Lower  steam  temperature.  Lower  steam  tempera¬ 
tures  could  also  result  in  significant  reductions  in  dis¬ 
tribution  system  heat  losses.  For  instance,  a  reduction 
in  supply  temperature  from  375° F  to  320° F  would  re¬ 
sult  in  more  than  15%  reduction  in  heat  losses.  The  re¬ 
duction  in  supply  temperature  could  be  made  up  of  a 
reduction  in  steam  superheat  as  well  as  reduced  steam 
saturation  temperature  due  to  reduced  pressure.  An 
analysis  similar  to  the  one  above  gives  the  present  value 
of  future  heat  savings  as  $1,433,013.  Since  this  might 
well  be  accomplished  with  little  or  no  investment,  it's 
very  attractive. 

CONCLUSIONS  AND  RECOMMENDATIONS 

Several  major  opportunities  exist  for  reducing  the 
heat  losses  from  the  Fort  Wainwright  heat  distribution 
system.  During  the  average  year  the  cost  of  heat  loss 
approaches  one  million  dollars.  The  two  most  promis¬ 
ing  alternatives,  adding  1  inch  of  insulation  to  the  steam 
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pipes  in  the  larger  utilidors  and  decreasing  the  steam 
temperature  from  375°  to  320°F,  were  found  to  offer 
the  potential  for  significant  savings. 

Each  of  these  alternatives  merits  further  study.  The 
results  of  this  work  are  limited  primarily  by  the  qual¬ 
ity  of  the  input  data.  More  detailed  studies  should  be 
performed  before  any  alternative  is  undertaken.  Thus, 
our  recommendations  will  deal  with  refinements  to 
this  study  which  would  yield  results  of  higher  confi¬ 
dence.  Our  recommendations  are  listed  below: 

1 .  Establish  more  accurate  estimates  as  to  the  thick¬ 
ness  and  condition  of  the  existing  pipe  insulation.  Lab¬ 
oratory  measurements  of  the  thermal  conductivity  of 
actual  samples  of  the  insulation  are  also  needed.  These 
tests  are  relatively  routine. 

2.  Determine  the  thermal  characteristics  of  the  soil 
system  around  typical  utilidors.  The  long-term  effects 
of  such  a  heat  source  within  the  soil  may  significantly 
effect  the  thermal  properties  of  the  soil  by  drying  it 
out. 

3.  Obtain  more  accurate  geometric  data  on  the 
utilidors.  These  data  would  include  exact  interior  and 
exterior  dimensions,  burial  depth  and  location  of  the 
various  utilities  within. 

4.  Make  temperature  measurements  within  utili¬ 
dors  to  estimate  the  accuracy  of  the  heat  loss  compu¬ 
tational  procedure.  Heat  flux  meters  could  also  be  in¬ 
stalled  in  representative  utilidors  for  further  confirma¬ 
tion. 

5.  Refine  the  methods  used  in  the  computational 
procedure  for  determining  the  heat  losses.  Use  several 
full  numerical  models  to  investigate  selected  utilidor 
configurations. 

6.  Investigate  the  consequences  of  lowering  steam 
pressure  (saturation  temperature)  and/or  steam  super¬ 
heat  to  obtain  lower  overall  steam  temperatures. 

7.  Make  preliminary  estimates  of  the  cost  of  imple¬ 
menting  the  suggestions  put  forth  in  this  study.  If 
warranted,  make  more  detailed  estimates  of  the  attrac¬ 
tive  alternatives. 
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APPENDIX  A.  HEAT  DISTRIBUTION  SYSTEM  DATA 


Utllldor  Supply  pipe  Return  pipe 

Combination  Dimensions  (ft)  diameter  diameter  Length 


Number 

Width 

Height 

(in) 

(In) 

(ft) 

i 

1.00 

1.00 

1.50 

1.00 

100.0 

2 

1.50 

1.50 

0.00 

1.00 

130.0 

3 

1.50 

2.00 

2.00 

0.00 

1AO.O 

4 

2.00 

2.00 

1.25 

1.00 

35.0 

5 

2.00 

2.00 

2.00 

2.00 

1280.0 

6 

2.00 

2.00 

2.50 

2.00 

195.0 

7 

2.00 

2.00 

3.00 

1.25 

100.0 

8 

2.00 

2.00 

3.00 

1.50 

150.0 

9 

2.00 

2.00 

A.  00 

2.00 

170.0 

10 

2.00 

2.00 

5.00 

1.50 

60.0 

11 

2.00 

2.00 

10.00 

6.00 

35.0 

12 

2.00 

1.50 

2.00 

1.00 

30.0 

13 

2.00 

2.50 

2. 00 

1.00 

325.0 

1A 

2.00 

2.50 

2.00 

2.00 

275.0 

15 

2.00 

2.50 

6.00 

2.50 

620.0 

16 

2.00 

3.00 

6.00 

2.50 

620.0 

17 

2.00 

A. 00 

2.00 

1.50 

120.0 

18 

2.00 

A. 00 

3.00 

2.00 

275.0 

19 

2.A2 

2.A2 

2.50 

1.00 

190.0 

20 

2.50 

2.00 

2.00 

2.00 

150.0 

21 

2.50 

2.00 

6.00 

2.50 

25.0 

22 

2.50 

2.00 

6.00 

0.00 

30.0 

23 

2.50 

2.50 

2.00 

2.00 

65.0 

2A 

2.50 

3.00 

6.00 

A. 00 

225.0 

25 

2.50 

3.00 

3.00 

A. 00 

190.0 

26 

2.50 

3.50 

A. 00 

2.00 

30.0 

27 

3.00 

2.00 

8.00 

2.00 

415.0 

28 

3.00 

2.50 

3.00 

1.50 

60.0 

29 

3.00 

2.75 

3.00 

1.00 

100.0 

30 

3.00 

2.75 

A. 00 

2.00 

510.0 

31 

3.00 

3.00 

1.00 

1.00 

70.0 

32 

3.00 

3.00 

1.50 

1.50 

50.0 

33 

3. 00 

3.00 

2.00 

1.50 

80.0 

3A 

3.00 

3.00 

2.00 

1.25 

55.0 

35 

3.00 

3.00 

2.00 

2.00 

550.0 

36 

3.00 

3.00 

2.50 

1.00 

75.0 

37 

3.00 

3.00 

2.50 

1.50 

230.0 

38 

3.00 

3.00 

3.00 

1.50 

635.0 

39 

3.00 

3.00 

3.00 

2.00 

575.0 

AO 

3.00 

3.00 

A.  00 

2.00 

395.0 

Al 

3.00 

3.00 

6.00 

2.00 

305.0 

A2 

3.00 

3.00 

6.00 

3.00 

670.0 

A3 

3.00 

3.00 

8.00 

2.00 

620.0 

AA 

3.00 

3.00 

8.00 

4.00 

270.0 

A5 

3.00 

3.00 

8.00 

A.  00 

210.0 

A6 

3.06 

0.00 

6.00 

0.00 

0.0* 

A7 

3.00 

3.50 

1.25 

1.25 

330.0 

A8 

3.00 

3.50 

1.50 

1.50 

60.0 

A9 

3.00 

3.50 

2.00 

1.00 

260.0 

50 

3.00 

3,50 

2.00 

1.50 

805.2 

51 

3.00 

3.50 

2.00 

2.00 

60.0 

52 

3.00 

3.50 

2.50 

1.50 

965.0 

53 

3.00 

3.50 

3.00 

1.25 

330.0 

5A 

3.00 

3.50 

3.00 

1.50 

935.0 

55 

3.00 

3.50 

3.00 

2.00 

675.0 

56 

3.00 

3.50 

3.50 

2.00 

75.0 

57 

3.00 

3.50 

A. 00 

2.00 

6360.0 

58 

3.00 

3.50 

5.00 

2.50 

A5.0 

59 

3.00 

3.50 

6.00 

2.00 

525.0 

60 

3.00 

3.50 

6.00 

3.00 

650.0 

61 

3.00 

3.50 

8.00 

A. 00 

700.0 

62 

3.00 

A. 00 

2.00 

1.00 

40. 0 

63 

3.00 

A. 00 

3.00 

2.00 

270.0 

6A 

3,00 

A. 00 

4.00 

3.00 

465.0 

65 

3.00 

A. 00 

6.00 

3.00 

50.0 

66 

3.00 

5.00 

2.00 

2.00 

65.0 

67 

3.50 

3.00 

2.50 

1.00 

160.0 

68 

3.50 

3.00 

4.00 

2.00 

150.0 

69 

3.50 

3.00 

6.00 

2.50 

150.0 

70 

3.50 

3.00 

6.00 

3.00 

520.0 
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Utilidor 

Combination  Dimensions  (ft) 


Number 

Width 

Height 

71 

3.50 

3.50 

72 

3.50 

3.50 

73 

3.50 

3.50 

74 

3.50 

4.50 

75 

4.00 

3.00 

76 

4.00 

3.00 

77 

4.00 

4.00 

78 

4.00 

4.00 

79 

4.00 

4.00 

80 

4.00 

4.00 

81 

4.00 

4.00 

82 

4.00 

5.00 

83 

4.00 

5.00 

84 

4.00 

5.00 

85 

4.00 

5.00 

86 

4.00 

5.00 

87 

4.00 

5.00 

88 

0.00 

0.00 

89 

4.00 

5.00 

90 

4.00 

5.00 

91 

4.00 

5.00 

92 

4.00 

5.00 

93 

4.00 

5.00 

94 

4.00 

5.00 

95 

4.00 

5.00 

96 

4.00 

5.00 

97 

4.00 

5.00 

98 

4.50 

4.00 

99 

4.50 

5.00 

100 

4.50 

5.00 

101 

4.50 

5.00 

102 

0.00 

0.00 

103 

4.50 

5.00 

104 

4.50 

5.00 

105 

4.50 

5.00 

106 

4.50 

4.50 

107 

4.50 

4.50 

108 

4.50 

4.50 

109 

4.50 

4.50 

110 

4.50 

4.50 

111 

5.00 

4.00 

112 

5.00 

4.75 

113 

5.00 

5.00 

114 

5.00 

5.00 

115 

0.00 

0.00 

116 

5.00 

5.00 

117 

5.00 

5.00 

118 

0.00 

0.00 

119 

5.00 

5.00 

120 

5.00 

5.00 

121 

5.00 

5.00 

122 

5.00 

5.00 

123 

5.00 

5.00 

124 

5.00 

5.00 

125 

5.00 

5.00 

126 

5.00 

5.00 

127 

0.00 

0.00 

128 

5.00 

5.00 

129 

5.00 

5.00 

130 

0.00 

0.00 

131 

5.00 

5.00 

132 

5.00 

5.00 

133 

5.00 

5.00 

134 

0.00 

0.00 

135 

5.00 

5.00 

136 

0.00 

0.00 

137 

0.00 

0.00 

138 

5.00 

5.00 

139 

5.00 

5.00 

Supply  pipe 
diameter 
(in) 

Return  pipe 
diameter 
(in) 

Length 

(ft) 

1.50 

1.50 

70.0 

6.00 

3.00 

420.0 

3.00 

4.00 

50.0 

6.00 

4.00 

520.0 

3.00 

2.00 

500.0 

3.50 

2.50 

345.0 

3.00 

2.00 

90.0 

4.00 

2.00 

515.0 

6.00 

2.00 

410.0 

6.00 

3.00 

510.0 

8.00 

2.00 

100.0 

1.00 

1.00 

120.0 

2.00 

0.00 

130.0 

0.00 

2.00 

290.0 

3.00 

1.50 

100.0 

4.00 

2.00 

260.0 

4.00 

2.50 

240.0 

0.00 

6.00 

0.0 

5.00 

3.50 

400.0 

6.00 

2.50 

40.0 

6.00 

3.00 

370.0 

6.00 

4.00 

270.0 

6.00 

0.00 

200.0 

8.00 

3.00 

730.0 

8.00 

4.00 

1015.0 

10.00 

4.00 

180.0 

10.00 

6.00 

1200.0 

6.00 

2.00 

190.0 

3.00 

1.00 

320.0 

3.00 

2.00 

190.0 

3.50 

2.00 

145.0 

0.00 

1.25 

0.0* 

5.00 

2.00 

230.0 

6.00 

3.00 

555.0 

10.00 

4.00 

365.0 

4.00 

2.00 

770.0 

6.00 

2.00 

240.0 

8.00 

4.00 

960.0 

10.00 

4.00 

170.0 

12.00 

6.00 

825.0 

8.00 

4.00 

460.0 

12.00 

8.00 

450.0 

1.00 

1.00 

130.0 

1.00 

1.00 

410.0 

8.00 

4.00 

0.0* 

1.25 

1.25 

425.0 

1.50 

1.25 

170.0 

6.00 

3.00 

0.0* 

2.00 

1.00 

400.0 

2.00 

2.00 

310.0 

3.00 

1.25 

190.0 

3.00 

1.50 

125.0 

3.00 

2.00 

300.0 

4.00 

1.50 

155.0 

4.00 

2.00 

2315.0 

4.00 

2.00 

400.0 

2.00 

0.00 

0.0* 

6.00 

2.50 

1115.0 

6.00 

3.25 

355.0 

2.00 

1.25 

0.0* 

6.00 

3.00 

500.0 

6.00 

4.00 

120.0 

2.50 

1.25 

2800.0 

8.00 

0.00 

0.0* 

8.00 

3.00 

410.0 

1.50 

1.50 

0.0* 

0.00 

1.25 

0.0* 

8.00 

4.00 

6285.0 

8.00 

4.00 

400.0 
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Combination 

Number 


Utilidor 
Dimensions  (ft) 
Width  Height 


Supply  pipe 
diameter 


Return  pipe 
diameter 


140 

0.00 

0. 

141 

0.00 

0. 

142 

5.00 

5. 

143 

5.00 

5. 

144 

5.00 

5. 

145 

0.00 

0. 

146 

5.00 

5. 

147 

0.00 

0. 

148 

5.00 

5. 

149 

5.00 

5. 

150 

5.00 

5. 

.51 

5.00 

5. 

1 52 

5.00 

5. 

53 

5.33 

6. 

.54 

0.00 

0. 

55 

5.50 

4. 

56 

5.50 

5. 

.57 

5.50 

5. 

.58 

5.50 

5. 

59 

5.50 

5. 

60 

5.50 

5. 

61 

6.00 

4J 

62 

6.00 

5.' 

63 

6.00 

6.i 

64 

0.00 

O.i 

65 

6.00 

6.1 

66 

6.00 

6.1 

67 

6.00 

6.i 

68 

6.00 

b. 

69 

6.00 

6. 

70 

6.00 

6.: 

71 

6.00 

7.1 

72 

6.00 

7.1 

73 

6.00 

7.1 

75 

6.0 

76 

0.0 

77 

6.0 

78 

6.0 

79 

0.0 

80 

6.0 

81 

6.0 

82 

6.0 

83 

0.0 

84 

6.0 

85 

0.0 

86 

0.0 

87 

6.0 

88 

6.0 

89 

6.0 

90 

0.0 

91 

6.0 

92 

0.0 

93 

6.0 

94 

0.0 

95 

6.0 

96 

6.0 

r 


Combination 

Number 

Utilidor 
Dimensions  (ft) 
Width  Height 

Supply  pipe 
diameter 
(in) 

Return  pipe 
diameter 
(in) 

Length 

(ft) 

211 

0.00 

0.00 

5.00 

0.00 

0.0* 

212 

7.00 

9.00 

14.00 

5.00 

50.0 

213 

0.00 

0.00 

5.00 

0.00 

0.0* 

214 

7.00 

9.00 

2.50 

1.25 

40.0 

215 

0.00 

0.00 

8.00 

0.00 

0.0* 

216 

7.00 

9.00 

10.00 

2.00 

100.0 

217 

7.00 

9.00 

10.00 

6.00 

290.0 

218 

0.00 

0.00 

10.00 

5.00 

0.0* 

219 

7.00 

9.00 

10.00 

6.00 

320.0 

220 

0.00 

0.00 

12.00 

5.00 

0.0* 

221 

0.00 

0.00 

3.00 

0.00 

0.0* 

222 

7.00 

9.00 

10.00 

6.00 

290.0 

223 

0.00 

0.00 

14.00 

5.00 

0.0* 

224 

0.00 

0.00 

3.00 

0.00 

0.0* 

225 

7.00 

9.00 

10.00 

5.00 

80.0 

226 

0.00 

0.00 

14.00 

6.00 

0.0* 

227 

7.00 

9.00 

12.00 

5.00 

90.0 

228 

7.00 

9.00 

12.00 

5.00 

280.0 

229 

0.00 

0.00 

5.00 

0.00 

0.0* 

230 

7.00 

9.00 

8.00 

1.50 

40.0 

231 

7.00 

9.00 

10.00 

5.00 

90.0 

232 

7.50 

6.00 

20.00 

10.00 

450.0 

*  Note:  In 

instances 

where  the 

utilidor  dimensions 

and  length 

are  zero, 

pipes  are  contained  in  preceeding  utilidor.  In  addition  3  pipes  were  in¬ 
advertantly  left  out  of  the  computer  input  file.  This,  unfortunately,  was  net 
noticed  until  all  calculations  were  complete.  Since  the  error  caused  by 
this  mistake  is  of  the  order  of  less  than  1/2%,  the  calculations  were  not 
redone. 
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APPENDIX  B:  UTILIDOR  HEAT  LOSS  PROGRAM 


INITIATE  CONSTANTS 
B0  =  4 

PI=3.14 159 
R I THK  = 1 
SITHKxl 
S0ILC=1.5 
T  CNOST  =  190 
THCC=G.8U 
TST£AM=37b 
TSURP=25. 7 
U  THK  =  •  5 

SITHK=S1THK/12. 

RITHK=RITHK/12. 

WRITE  INTROOUCT ION  ANO  HEADINGS 
WRITE(6.b> 

5  F0RMAT(1X,*DATA  OU7°UT*/// 

45X , 'THIS  PROGRAM  CALCULATES  THE  HEAT  LOSS  PROM* 

4/»  UTITLIOOR  SECTIONS  WITH  THE  VAR  I ALbES  bELOW.* 

i//*  independent:*. ri9,*x  =utiliddr  height  <ft>* 

4/18X,*T  -UTILIDOR  WIDTH  (FT)* 

1/ 18  X. *  L  =UTILinoR  LENGTH  (FT)* 

4/18X, *SPD  ^SUPPLY  PIPE  DIAMETER  (1N>* 

»/18X.*RPD  ^RETURN  PIPE  DIAMETER  (IN)* 

4/18X«*TCnDT=RETURN  CONDENSATE  TEMPERATURE* 
i/18X«*TSTEM-SUPPLY  STEAM  TEMPERATURE* 

4/18X.*TSURF=GR0UND  SURFACE  TEMPERATURE* 

4//*  DEPENDENT:*, T la, • TBULK=bULK  AIR  TEMPERATURE* 

4/ 18X  « *TCONC-INTER IOR  CONCRETE  TEMPERATURE* 

4/18X,*TC0NX=EXTERI0R  CONCRETE  TEMPERATURE* 

4/18X.*TINSS=SUPPLY  INSULATION  TtMPERATURE* 

4/18X,*TINSR=RETURN  INSULATION  TEMPERATURE*//) 

WRITE (fa. b) 

6  FORMAT ( 1H  /3X.*X*,BX.*Y*,4X,*SPD*.3X,*RPD'.BX,11HTSURF  T  ST  EM, 

SIX  .S2HTCNDT  TINSS  TINSR  TfcULK  TCONC  TCONX  HEAT  LOSS, 

44X,*L*,6X,*HEAT  LOSS • / 2 X . 3 3H ( f T )  (FT)  (IN)  (IN)  (F), 

4  3X , 4  7H ( F )  (F>  (F)  (F)  (F)  (F)  (F)  , 

4*(BTU/HR*FT)* ,?X, *(FT) • , 5 X , • ( B TU/H R > */lX»116(*-*>) 

READ  IN  NJMBER  OF  SECTIONS 
REAO(S,7)  N 

7  FORMAT  (  14  ) 

DO  1000  1=1. N 

READ  IN  UTILIDOR  AND  PIPE  DIMENSIONS 
R  EAD ( 5.8 )  X#Y  ,SPD,RPD,XL 

8  FORMAT (IX,F4.2,3F6.t,F7.1) 

CHECK  FOR  INPUT  ERROR 
CHECK=X»Y*XL 

I F ( CHECK )  10,9,11 

9  IF(X.EQ.O.AND.Y.EQ.O«AND.XL.EQ.C>  GO  TO  11 

10  WRITr(6,ll)  I ,X,Y,S?D,RPO,XL 
GOTO  1000 

11  FORMATUHO,  'UTILIDOR  DIMENSION  INPUT  ERROR  ON  LINE  *,I4 
4,*.  DATA  ENTERE0=*.4F5.2,F7.1> 

DETERMINE  TYPE  OF  SECTION 
IF(RPO.EQ.O.ANO.SPO.GT.O)  GOTO  15 
IF(SPD.EQ.O.ANO.RPO.GT.O)  GOTO  1 A 
IF(SPD.GT.O.ANO.RPO.GT.O)  GOTO  13 
WRITE(6,12>  I,X,Y,SPD*RPO»XL 
GOTO  1000 

12  FORMATC1X, ‘ERROR  IN  PIPE  DIMENSION  DATA  ON  LINE  NUMBER  • ,! A 
4,«.  DATA  LINE=*,1X,FA.2,3F6.2,F7. 1) 

13  KASE=1 
GOTO  16 

1A  KASE=0 
GOTO  16 

15  KASE=-1 

16  CONTINUE 

CHECK  FOR  MULTIPLE  PIPES  ANO  SAVE  UTILIDOR  DIMENSIONS 
I F ( X )  21,21,22 
21  X  =  S  AVE  X 
Y  =  SAVE  T 
XL=SAVEL 

UTILIDOR  HEAT  LOSS  PROGRAM 

THIS  PROGRAM  CALCULATES  THE  HEAT  LOSS  FROM  PIPES  WITHIN 
BURIED  CONCRETE  UTIL1DORS.  ENTER  UTILIDOR  DATA  INTO  DATA¬ 
FILE  "DATAUT".  THE  FIRST  LINE  MUST  PE  THE  TOTAL  NUMBER  OF 
UTILIDOR  SECTIONS  TO  BE  EVALUATED  (14  FORMAT).  EACH  SUCCESSIVE 
ENTRY  LINE  MUST  CONTAIN  ,  IN  ORDER,  THE  UTILIDOR  HEIGHT  4 
WIDTH,  THE  STEAM  (SJPPLY)  PIPE  DIAMETER,  THE  CONDENSATE  (RETURN) 
PIPE  DIAMETER  ANO  FINALLY  THE  UTILIDOR  LENGTH  (ALL  PFR  FORMAT  »(’), 
THERE  MAY  BE  ANY  NUMBER  OF  PIPES.  THE  FOLLOWING  SYSTEM  IS 
EMPLOYED  FOR  ENTERING  THE  DIFFERENT  QUANTITY  AND  TYPE  OF  PIPES! 

1  PIPE:  ENTER  UTILIDOR  DIMENSIONS  AND  DIAMETER  OF  PIPF  IN 

CORRECT  CATEGORY  (SUPPLY  OR  RETURN).  ENTER  ZERO  (0) 
FOR  THE  MISSING  PIPE  DIAMETER. 

2  PIPES!  ENTER  UT1LIOOR  DIMENSIONS  AND  DIAMETERS  OF  PIPES 

IN  CORRECT  ORDER  (SUPPLY  THFN  RETURN). 
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MORL  THAN 

2  PIPES!  ENTER  MRST  T.O  (?)  PIPES  AS  ABOVF.  THEN  ON  A 
NEW  DATA  LINE  ENTER  ZEROS  (0)  FOR  ALL  UTIL1D0R 
DIMENSIONS  AND  ENTER  DIAMETER  OF  PIPE  IN  CORRECT 
CATCGCRY « 

EERO  DIMENSIONS  FOR  UTILIDDR  CAN  3F  AOPf C  1 NnFF I NF  TL  Y  TO 
INDICATE  ANY  NJ^bER  OF  PIPES. 

FOLLOWING  IS  A  LIST  OF  THE  SIGNIFICANT  PROGRAM  VARIABLES. 

60=  bJKlAL  DEPTH  OF  UTIL1D0R  (FT) 

EFAIRC=  EFECIIVE  AIR  CONDUCTIVITY  (u TU /HR *F T *F ) 

EFUD=  EFFECTIVE  UTILIDCK  DIAMETER  (FT) 

ORSHOFr  C-RASHOF  NUMBER 

HL  0  S  S  =  HEAT  LOSS  (STU/HK*FT) 

HTCA:  HEAT  TRANSFFR  COEF1CIENT  OF  AIR  ( BTU/hR * F T • F ) 

“TCI=  HEAT  TRANSFER  COEFICIENT  OF  INSULATION  ( RTU/HR  *F  T  *F ) 

PC  KU=  PERIMETER  OF  UTILIDOR  (FT) 

PR  ANT  L=  PRANDTL  NUMBER 

RAIRr  THERMAL  RESISTANCE  OF  AIR  (HR*FT*F/HTU) 

RC0NC=  THERMAL  RESISTANCE  OF  CONCRETE  ( HR *F T *F /P TU ) 

RINS(S)(R)=  THERMAL  RESISTANCE  OF  INSJLATION  ( SUPPLY )< R E T URN ) 
(HR  *F  T  *F/3  TU  > 

KlRr  RETURN  PIPE  INSULATION  RADIUS  (FT) 

R1THK=  RETURN  PIPE  INSULATION  THICKNESS  (FT) 

RPU=  RETURN  PIPE  DIAMETER  (IN) 

RPR=  .RETURN  PIPE  RADIJS  (FT) 

R  SO 1L  =  THERMAL  RESISTANCE  OF  SOIL  (HR  *  FT *F /HTu > 

BELT  A  =  EFFECTIVE  P  IPE./U  T  ILIOOR  DISTANCF  (FT) 

S I R  =  SUPPLY  PIPE  RADIUS  (FT) 

S I T HK  =  SUPPLY  PIPE  INSULATION  THICKNESS  (IN) 

S  0  I LC=  SOIL  CONDUCTIVITY  ( BT U /HR . F T»F > 

SPO=  SUPPLY  PIPE  DIAMETER  (IN) 

SPR=  SUPPLY  PIPE  RADIUS  (FT) 

T  BULK  =  BULK  AIR  TEMPERATURE  (F> 

T  CNDST  =  CONDENSATE  TEMPERATURE  (F> 

TC0NC=  INICRIOR  CONCRETE  TEMPERATURE  (F) 

THCAIR:  THERMAL  CONDUCTIVITY  OF  AIR  ( B  T  U /HR  *F  T  *F ) 

THCC=  THERMAL  CONDUCTIVITY  OF  CONCRETE  ( BT U /HR *F T *F ) 

T HCS I =  THERMAL  CONDUCTIVITY  OF  SUPPLY  INSULATION  ( BTU/HR* F T »F > 
THCRI=  THERMAL  CONDUCTIVITY  OF  RETURN  INSULATION  ( PT U /HR* F T *F ) 
TINS(S) (R)r  INSULATION  TEMPERATURE  ( SUPPL Y ) ( RE T UR N )  CF) 

TSTL  AM=  STEAM  TEMPERATURE  (F) 

T S JR F  =  GROUNO  SURFACE  TEMPERATURE  «F> 

TX=  EXTEkIOR  CONCRETE  TFMPERATURF  (F) 

JTHK=  UTILIOOR  THICKNESS  (FT) 

V  A I R  =  DYNAMIC  VISCOSITY  OF  AIR  (FT**?/SEC) 

X:  UTILIDDR  HEIGHT 

XL=  UTILIOOR  LENGTH  (FT) 

Y=  JT ILIOOR  WIDTH  (FT) 

DIMENSION  0(000), QXL(OOO) 

CALL  C3NTRL(3,*UATAJT«,b> 

CALL  CONTRL  (3  ,  »  TO  TAL  (.  •  ,  6  ) 
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MULT  =  MUL 
GO  TO  25 
SAVEXrX 
SAVEYrY 
SAVEL=XL 
MULT  =  0 
S  I  R  =  0 
R  I R  =  0 
HLOSS  X  =  0 
LASTrC 


T*J 


C 

c 


C 

r 


CONVERT  UNITS  AND  DIAMETERS  TO  RADII 

SPD=SPD/12 

RPD-RPD/12 

SPR=SPD/2 

RPR=RPD/2 

IF(KASE)  27,28,2b 

26  R  IR  =RPR  +R  I  T  HK 

27  SIR=SPR*SITHK 
GO  TO  29 

28  RIR-RPR*RITHK 

29  RIO:RIR*2. 

SI0=S1R*2. 

GUESS  AT  CONCRETE  AND  INSULATION  TEMPERATURES 

TINSS:TSTEAM-10D 

TINSR=TCNDST-100 

TC0NC=55 

ITNUMrE 

IF(KASE)  33,00,55 
33  T INS=T INSS 
GOTO  60 
00  TINS=TINSR 
GOTO  60 

55  UNS=<T INSS*SIU*TINSK*RIOT/<SID*RIO) 

60  Tt)ULK  =  (TINS*TC0NC>/2. 

2  TNUM  =  1 TNJM*I 


CALCULATE  THE  THERMAL  RESISTANCE  FOR  EACH  MAT*  L 
?HCSI=.0221*O.13E  -5*<TSTEAM*TINSS)/2 
THCRI=.0221*0.13E  -5* <TCNOST*T INSR )/2 
IF(KASE)  63*60,63 

63  RINSS=AL0GTSIR/SPR)/(2*PI*THCSI ) 

T F ( K AST  >  65,65,60 
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64  RINSR=ALQGCKIk/RPR)/C2.*PI  *  I H  CR 1  ) 

65  PER'J=2,  » ( X»Y) 

RCONC-UTHK/<PEkU*ThCC) 

RS01L  =  l/CS0ILC*l.(.Ho*<AL  9G  1  3  <  :-')/x»  1  .  )>••«-.  *'>)•!<  KJ/YWC-.  0  76  )  )  ) 

PRANTL=C.71dM9-TBUL<*1.275f  -4 

EFUD=PFRU/°I 

THCAIR=G.01319*T6ULK‘2.5C  -5 
V  AI  R  =  T  oULK  *5. 4E  -7»1.2P24f  -A 
DLLTAr<LFUD-CSID*RI J> )/2. 

GRSH0F  =  32.2*ABSCT  INS-TCONO  .(DELTA  ••«.  >/<<  VA  1R*,2.).<TPULK»450.7>  ) 

EF  AIRC=ThCA1R*0.4*(GRShGF  **‘R  A  NT  L  )  »  *o*2 
IF (<ASE)6b,67  *70 

66  RAIR=AL0G<EFUD/SID)/C2.*PI*LFAIRC) 

RINS-RIMSS 

GOTO  60 

67  RAlR=AL0uCEFUD/RI0>/<2.»PI»C''AIAC> 

R I NS=R I NSR 

GOTO  80 

70  HTC1=1./RIN5SMTCNPST-TINSR)/C<TSTEAM-TINSS)‘RINSR) 

HTCA=2*PI*EFAIRC*(1/ALQG(EFU0/SID)‘<TIN5R-TC0NC)/(CTINSS-TC0NC> 
&*ALOG<EFUO/RID>  > ) 

R1NS=1./HTC1 

RAIR=1*/HTCA 

SUM  THERMAL  RESISTANCES  AND  FIND  HEAT  LOSS. 

80  HTCT=1./CR1NS+RAIR+RC0NC+RS0IL> 

IFCKASE)  b7,66,R7 

86  MLOSS=HTCT*(TCNDST-TSURF) 

TINSR-TCNOST-CHLOSS»RINS) 

T I NS=T I NSR 

GOTO  90 

87  HLOSS  =  HTCT* < T STE AM-T SURF  ) 

TINSS=TSTEAF-(HLOSS*RINS) 

T I  NS  =  T I NSS 

90  TX=TSJRF,(HLOSS*RSOIL) 

TCONC=TX*CHLOSS*RCONC) 

IF(KASC)  94,94,93 

93  RAIRS=ALO&CEFUO/SIO)/(2.*PI*EFAIRC> 
RAIRR=ALOG(EFUD/RIO)/(2.*Pl*EFAIRC) 

TINSS  =  TSTEAM/(l»RlNSS/RAIRS)*TCONC/(  1  +  RA1KS/R1NSS) 
TINSR=TCNDST/(1*RINSR/RAIRR)*TC0NC/I1*RAIRR/RINSR> 


95 

96 

100 

101 

105 

no 

120 


DETERMINE  if  heat  loss  value  has  converged  and  continue 

ITERATIONS  IF  NECESSARY. 

DELTAH=AbS(HL0SS-HL3SSX) /HLOSS 

PERCNT-HLOSS*O.OQOQ1 

IF(OELTAH.LE.PERCNT)  GOTO  100 

HLCSSX-HLOSS 

IF<ITNUM.EQ.121  GOTO  95 

IFCKASE)  60,60,55 

URITEt6,9b  > 

FORMAT  (  5X, *AFTER  12  ITERATIONS  CQNVEKGIANCE  NOT  ATTAINED*/ 
12X, ‘ALTER  INITIAL  TEMPERATURES  IN  LINES  ABOVE*) 

GOTO  1000 
0  < I ) =HLOSS 

IF(LAST)  101,101,105 
LAST=1.0 

IFCKASE)  60,60,55 

FIND  HEAT  LOSS  OVER  PIPE  LENGTH 

OXLCI)=Q(I)*XL 

I F ( MULT )  120,120,110 

X  =  0 

Y  =  0 

XL  =  0 

S  =  SP0  *1 2 • 

R=RP0*12. 


C 


I  ITERATIONS  COMPLETE,  WRITE  RESULTS 

WRITE! 6, 140)  X , Y, S . R , T SUR F , T S TE AM , TCNDST , T I NSS , T I NSR .TPULK 
&TCONC,TX*S( I)  ,  XL  »  Q  XL ( I ) 

140  F0RMAT(1X,4(F5.2,1X),5X,4«F5.1,IX),5X,4CF5.1,1X),5X,F6.2, 
&2X,F6.1«2X,E12.5) 

180  XI=I/4. 

ISPACErlFIX(XI) 

IFUSPACE.EG.XI  )  GOTO  190 
GOTO  196 
190  VRITE(6,195> 

195  FORMAT  C 1H  ) 

196  XI=A6SCCI-40)/57.) 

I  TI  TLE  =  I  F  I  X  (  X  I  > 

IFCIT1TLE.EQ.XI)  WR1TEC6.6) 

1000  CONTINUE 


1200 

C 


1300 


SECTIONS  COMPLETE,  SUM  HEAT  LOSS 
SUM=0 . 

DO  1200  J  =  1 ,  N 
SUM-SUM*uXL< J) 

CONTINUE 

WRITE  OUT  TOTAL  HEAT  LOSS 
SUMYR=SUM*0. 00876 
WRITEC6.1300)  SUM , SUM  YR 

FORMAT  C 1  HO//, •  TOTAL  HEAT  LOSS  FROM  UTILIDDRS-**E12.4«* 
t/T33,E12.4, *  MbTU/YR  •  ) 

CALL  C0NTRLC4, ‘TOTALS*, 6) 

CALL  C0NTRL<4,‘UATAJT‘,5) 

CALL  EXIT 
ENO 


8 T U/HR  * 
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APPENDIX  C:  SAMPLE  INPUT  AND  OUTPUT  DATA  FILES 


OAT*  INPUT  FILE 

THIS  FILE  CONTAINS  THE  VARIABLES  REQUIRED  OF  THI  THl  MEAT  LOSS 
PROGRAM:  "GLOSS".  THE  FIRST  LINE  IS  THE  TOTAL  NUMBER  OF  SECTIONS 
TO  BE  EVALUATED.  EACH  COLUMN •  THEREAFTER.  CONTAINS  VALUES  FOR  THE 
VARIABLES  BELOV  R ESP E C T I VE LT : 

X  =UT1L1D0R  HEIGHT  ( FT  > 

Y  =UT IL l DOR  VIDTH  !FT> 

SPD  -SUPPLY  PIPE  DIAMETER  UN) 

SKO  ^RETURN  PIPE  DIAMETER  (IN) 

L  =UT I L I  DOR  LENGTH  (FT) 
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note:  all  comments  must  he  deleted  before  using 

THIS  FILE. 


X 

Y 

SPD 

SRD 

L 

1.0C 

1.00 

1.50 

1.00 

100.0 

1.50 

1.50 

0.00 

1.00 

130.0 

1  .50 

2.00 

2.00 

0.00 

140.0 

2.00 

2.00 

1  .25 

1.00 

35.0 

2.00 

2.00 

2.00 

2.00 

1280.0 

2.00 

2.00 

2.50 

2.00 

195.0 

2.00 

2.00 

3 . 0  G 

1.25 

100.0 

2.00 

2. 00 

3.00 

1.50 

150.0 

2.00 

2.00 

4.00 

2.00 

170.0 

2.00 

2.00 

5.00 
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